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Abstract. We define and study distributions in M'' that we call g— Normal. 
For q = I they are really multidimensional Normal, for q € (—1, 1) they have 
densities, compact support and many properties that resemble properties of 
ordinary multidimensional Normal distribution. We also consider some gen- 
eralizations of these distributions and indicate close relationship of these dis- 
tributions to Askey- Wilson weight function i.e. weight with respect to which 
Askcy- Wilson polynomials are orthogonal and prove some properties of this 
weight function. In particular we prove a generalization of Poisson-Mehler 
expansion formula. 



1. Introduction 

The aim of tliis paper is to define, analyze and possibly 'accustom' new distri- 
butions in R*^. They are defined with a help of two one-dimensional distributions 
that first appeared recently, partially in noncommutative context and are defined 
through infinite products. That is why it is difficult to analyze them straightfor- 
wardly using ordinary calculus. One has to refer to some extent to notations and 
results of so called g— series theory. 

However the distributions we are going to define and examine have purely com- 
mutative, classical probabilistic meaning. They appeared first in an excellent paper 
of Bozejko et al. |4] as a by product of analysis of some non-commutative model. 
Later they also appeared in purely classical context of so called one-dimensional 
random fields first analyzed by W. Bryc at al. in [Tj and [3l . From these papers we 
can deduce much information on these distributions. In particular we are able to 
indicate sets of polynomials that are orthogonal with respect to measures defined 
by these distributions. Those are so called q— Hermitc and Al-Salam-Chihara poly- 
nomials - a generalizations of well known sets of polynomials. Thus in particular 
we know all moments of the discussed one-dimensional distributions. 

What is interesting about distributions discussed in this paper is that many of 
their properties resemble similar properties of normal distribution. As stated in the 
title we consider three families of distributions, however properties of one, called 
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multidimensional g— Normal, are main subject of the paper. The properties of the 
remaining two are in fact only sketched. 

All distributions considered in this paper have densities. The distributions in 
this paper are parametrized by several parameters. One of this parameters, called 
q, belongs to (—1, 1] and for q = 1 the distributions considered in this paper become 
ordinary normal. Two out of three families of distributions defined in this paper 
have the property that all their marginals belong to the same class as the joint, 
hence one of the important properties of normal distribution. Conditional distri- 
butions considered in this paper have the property that conditional expectation of 
a polynomial is also a polynomial of the same order - one of the basic properties 
of normal distributions. Distributions considered in this paper satisfy Gebelein in- 
equality -property discovered first in the normal distribution context. Furthermore 
as in the normal case lack of correlation between components of a random vectors 
considered in the paper lead to independence of these components. Finally condi- 
tional distribution fc {x\y, z) considered in this paper can be expanded in series of 
the form fc {x\y, z) = Jm {x) YhLo hi{x)9i iVi ^) where /m is a marginal density, 
{hi} are orthogonal polynomials of Jm and gi (y, z) arc also polynomials. In par- 
ticular if fc {x\y,z) = fc {x\q) that is when instead of conditional distribution of 
X\Y, Z we consider only distribution of X\Y then gi (y) = hi (y). In this case such 
expansion formula it is a so called Poisson-Mehler formula, a generalization of a for- 
mula with hi being ordinary Hermite polynomials and /m {x) = exp(— x^/2)/-\/27r 
that appeared first in the normal distribution context. 

On the other hand one of the conditional distributions that can be obtained 
with the help of distributions considered in this paper is in fact a re-scaled and 
normalized (that is multiplied by a constant so its integral is equal to 1) Askey- 
Wilson weight function. Hence we are able to prove some properties of this Askey- 
Wilson density. In particular we will obtain a generalization of Poisson-Mehler 
expansion formula for this density. 

To define briefly and swiftly these one-dimensional distributions that will be 
later used to construct multidimensional generalizations of normal distributions, let 
us define the following sets 



Let us set also m + S (q) = {x = m + y,y ^ S (q)} and m + S ((j) = (mi + S (q)) x 
. . . X (md + S {q)) if m = (mi, . . . , md). Sometimes to simplify notation we will use 
so called indicator functions 



The two one-dimensional distributions (in fact families of distributions) are given 
by their densities. 

The first one has density: 




[-2/yr^,2/^/r^ if \q\<i 

{-1,1} if q = -l 




X e A 
x ^ A 



(1.1) 




/jv (a;|^) = 



277^4 - (1 - q)x- 



f n ((1 + - (1 - 11(1 - q'+')Isi,) {x) 
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defined for |g| < 1, x e M. We will set also 

(1.2) /jv(x|l) = -^exp(-xV2). 



For q = — I considered distribution does not have density, is discrete with two equal 
mass points at S (—1). Since this case leads to non-continuous distributions we will 
not analyze it in the sequel. 

The fact that such definition is reasonable i.e. that distribution defined by 
Jn {x\q) tends to normal N (0, 1) as q — > 1~ will be justified in the sequel. The 
distribution defined by /jv {x\q) , —1 < q < 1 will be referred to as g— Normal 
distribution. 

The second distribution has density: 



(l-3a) fcN{x\y,p,q) = ^ 



27r^4 - (1 - 5)a;2 
(1.3b) 

- (i-pV)(i-gfe+i)((i + gfe)^-(i-g)xV) 

M " ~ " + P^q^^)xy + (1 - q)pHx^ + y^)q^'^ ^'"^ 

defined for \q\ < l,\p\ < 1, x £ R, y £ S (q). It will be referred to as {y, p, g)— Conditional 
Normal, distribution. For g = 1 we set 

f r I i\ 1 I i^- pyf\ 

fcN {x\y, p, 1) = -j===== exp 



v/2^(l-p2) \^ 2(l-p2)y 

(in the sequel we will justify this fact). Notice that we have fcN {x\y,0,q) = 
Jn {x\q) for aU y e S (q). 

The simplest example of multidimensional density that can be constructed from 
these two distribution is two dimensional density 

9 {x, y\p, q) = fcN {x\y, p, q) /n {y\q) , 

that will be referred to in the sequel as N2 (0,0, 1, l,p\q). Below we give some ex- 
amples of plots of these densities. One can see from these pictures how large and 
versatile family of distributions is this family 
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P = --6, q = -.7 

It has compact support equal to S (q) x S (q) and two parameters. One playing 
similar role to parameter p in two-dimensional Normal distribution. The other 
parameter q has a different role. In particular it is responsible for modality of the 
distribution and of course it defines its support. 

As stated above, distribution defined by /^v {x\q) appeared in 1997 in 4 in 
basically non-commutative context. It turns out to be important both for classical 
and noncommutative probabilists as well as for physicists. This distribution has 
been 'accustomed' i.e. equivalent form of the density and methods of simulation of 
i.i.d. sequences drawn from it are e.g. presented in '18|. Distribution fcN, although 
known earlier in nonprobabilistic context, appeared (as an important probability 
distribution) in the paper of W. Bryc [l] in a classical context as a conditional dis- 
tribution of certain Markov sequence. In the following section we will briefly recall 
basic properties of these distributions as well as of so called g— Hermite polynomials 
(a generalization of ordinary Hermite polynomials). To do this we have to refer to 
notation and some of the results of series theory. 

The paper is organized as follows. In section 2 after recall some of the results of 
g— series theory we present definition of multivariate g— Normal distribution. The 
following section presents main result. The last section contains lengthy proofs of 
the results from previous section. 

2. Definition of multivariate q-Normal and some related distributions 

2.1. Auxiliary results. We will use traditional notation of q— series theory 
i.e. [0], = 0; [n]^ = l + q + ...+ q-' = [n]^\ = iTU ^ with [0]^! = 1, E] ^ 

n > k > 

[n-k]^\[k]^\ ' — — _ \yQ useful to use so called q— Pochhammer 

, otherwise 

symbol for n > 1 : (a|q)„ = n"=o^ (l " > with (a|g)o 1 , (ai, 02, . . . , afc|q)„ = 
riiLi {'^i\q)n- Often (a|g)„ as well as (ai,a2, . . . ,ak\q)n will be abbreviated to (a)„ 
and (ai, 02, ... , ak)^ j if it will not cause misunderstanding. 
It is easy to notice that (?)„ = (1 — 9)" [?^]^! and that 

* I , otherwise 
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Let us also introduce two functionals defined on functions g : R — > C, 

"2 I 2 '2. I '2, 

UWl^ / 1.9 (2^)1 fN{x)dx, \\g\\cL^ / IfWI fcN {x\y,p,q)dx 



'L ~ J llyllCL 

and sets: 

Liq) ^ {,g:K^C: ||5||^ <oo}, 

CL{y,p,q) - {g:R^C:\\g\\cL<^}- 

Spaces {L and {CL (y, p, q) , || are Hilbert spaces witli tlie usual def- 

inition of scalar product. 

Let us also define the following two sets of polynomials: 

-the g— Hermite polynomials defined by 

(2.1) Hn+i{x\q) = xHn{x\q) - [n]qHn-i{x\q), 

for n>l with H^i{x\q) = 0, Ho{x\q) = 1, and 

-the so called Al-Salam-Chihara polynomials defined by the relationship for 
n>0: 

(2.2) Pn+i{x\y,p,q) ^ {x ~ pyq''-)Pn{x\y, p,q) - (1 - p^q''^'^)[n]qPn-i{x\y, p,q), 

with P_i {x\y, p, q) = 0, Pq {x\y, p, q) = 1. 

Polynomials (|2.ip satisfy the following very useful identity originally formulated 
for so called continuous g— Hermite polynomials hn (can be found in e.g. Thm. 
13.1.5) and here below presented for polynomials H„ using the relationship 

(2.3) /i„ {x\q) = (1 - q)"/' i/„ (7r=l«) ' ^ 1' 

inin(n,m) ^ ^ j. ^ 

(2.4) H,AMq)Hrn{x\q)= ^ " [j]^lHn+^^2k{x\q). 

— 7 7 

It is known (see e.g. yij) that g— Hermite polynomials constitute an orthogonal 
base of L (g) while from [3] one can deduce that {P„ (a^ly, p, g)}„>„2 constitute an 
orthogonal base of CL (y, p, g). Thus in particular = ^s(q)^^ i^lyi P^q) fcN {x\y, p,q) 
= E {X\Y — y) — py. Consequently, if F has also g— Normal distribution, then EXY 
= P- 

It is known (see e.g. [7] formula 13.1.10) that 

(2.5) |i?„(x|g)| < W„(g)(l-g)-"/^ 
where 



(2.6) (g) = 

i=0 

We will also use Chebyshev polynomials of the second kind t/„ (x), that is 
Un (cos 6) — ^"^l^^^g"^^ and ordinary (probabilistic) Hermite polynomials Hn {x) i.e. 
polynomials orthogonal with respect to exp(— a;^/2). They satisfy 3— term re- 
currences: 

(2.7) 2xUn{x) = Un+l (x) + Un-l (x) , 

(2.8) xHn{x) = Hn+1 (x) + nHn-1 
with U-i (x) = H_i{x) = 0, Uo (x) = Hi (x) = 1. 
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Some immediate observations concerning g-Normal and (y, p, g)— Conditional 
Normal distributions are collected in the following Proposition: 

Proposition 1. 1. fcN{x\y,0,q) = fN{x\q). 

2. Vn > : H„ {x\0) = t/„ (x/2) , H„ (x|l) = i/„ (x) . 

3. Vn > : Pn{x\y,0,q) = Hn{x\q), Pn{x\y,p,l) = {l-p^r^^H„ 

P„ {x\y, p, 0) - Un {x/2) - pyU„^i {x/2) + p'U^^2 {xl2) . 

4. fN{x\Q) = ^V4- a;2/<_2,2> {x) , (x|q) — > -^exp(-x^/2) point- 

wise. 

5. fcN{x\y,p,0) - 2.((i-p^)^-Ki+l'S^P-(^-+.-)) ^<-^-^> ' fcN{x\y,p,q) 

Proof. 1. Is obvious. 2. Follows observation that (|2.ip simplifies to (|2.7p and 
(j2.8p for (J = and q = 1 respectively. 3. First two assertions follow either direct 
observation in case of P„ {x\y, p, 0) or comparison of (|2.2I) and (|2.8p considered for 
X — (x — py)/ y/l — p2 and then multiplication of both sides by (l — p2j("+i)/2 ^ 
Third assertion follows following observations: P_i {x\y, p, 0) = 0, Pq {x\y, p, 0) = 1, 
Pi {x\y, p,0) = X - py , P2 {x\y, p, 0) = x{x - py) - (l - p^) , Pn+i {x\y, p, 0) 
xPn {x\y, p, 0) — {x; y, p, 0) for n > 1 which is an equation p.7p with x replaced 
by x/2. 

4. 5. First assertions are obvious. Rigorous prove of pointwise convergence of 
respective densities can be found in work of [9] . To support intuition we will sketch 
the proof of convergence in distribution of respective distributions. To do this 
we apply 2. and 3. and see that \/n > 1 (x\q) — > i?„ (x) , and Pn(x\y,p,q) 



(1 - p^r'^H^ 



as q — !■ 1 . Now keeping in mind that families 

(a;|g)}„>g and {-Pn (xjy, p, q)}->Q are orthogonal with respect to distributions 
defined by respectively /at and fc n we deduce that distributions defined by Jm and 
fcN tend to normal -/V (0, 1) and A^(p?/,(l — p^)) distributions weakly as q — > 1~ 
since both N (0, 1) and N (^py, (l — p^)) are defined by their moments, which are 
defined by polynomials iJ„, and P„. □ 



2.2. Multidimensional g— Normal and related distributions. Before we 
present definition of the multidimensional g— Normal and related distributions, let 
us generalize the two discussed above one-dimensional distributions by introducing 
(m, CT^, g)— Normal distribution as the distribution with the density /nHx — m) /a\q)/a 
for TO e M, CT > 0, g e (-1, 1]. That is if X (to, ct^, g) -Normal then {X - m)/a 
^ g— Normal. 

Similarly let us extend definition of (y, p, g)— Conditional Normal by introduc- 
ing for TO e M, cr > 0, g e (—1,1], |p| < 1, (m, cr^, y, p, g)-Conditional Normal 
distribution as the distribution whose density is equal to fcN {{x — m)/a\y, p, g) /a. 

Let m, cr e M'' and p S (-1, 1^'^ , g S (-1, 1]. Now we are ready to introduce 
a multidimensional g— Normal distribution Nd (m, cr^, p|g) . 
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Definition 1. Multidimensional q— Normal distribution Nd (m, cr"^ , p\q) , is the 

continuous distribution in M."^ that has density equal to 

9 (^x|m, cr = /^r — \q j || fcN ^ 1 — ,P^,q j 

where x = {xi,. . .,Xdf , m = (mi, . . . , m^) , (t^= (aj, . . . , cr^) , p = (p^, . . -^Pd-i)- 

As an immediate consequence of the definition we see that supp(-/Vd(m, <T^|g)) 
= m + S (g) . One can also easily see that m is a shift parameter and cr is a scale 
parameter. Hence in particular EX = m. In the sequel we will be mostly concerned 
with distributions Nd{0, 1, p\q). 

Remark 1. Following assertion 1. of Proposition{J\we see that distribution Nd (0, 1, 0\q) 
is the product distribution of d i.i.d. q— Normal distributions. Another words "lack 
of correlation means independence" in the case of multidimensional q~ Normal 
distributions. More generally if the sequence p= {pi, . . . , Pd^i) contain, say, r 
zeros at, say, positions ti, . . . ,tr then the distribution of Nd {0,1, p) is a prod- 
uct distribution of r + 1 independent multidimensional q— Normal distributions: 
Nt, {0,l,{p„...,p,^_,)) ,...,Nd-t,. {0,l,iPt^+„...,pJ). 

Thus in the sequel all considered vectors p will be assumed to contain only 
nonzero elements. 

Let us introduce the following functions (generating functions of the families of 
polynomials) : 

(2.9) ^ix,t\q) = J2TT-^H,{x\q), 

i=o ^^J?- 

oo 

(2.10) Tix,t\y,p,q) = ^ P,(a;|y,p,g). 

■1=0 

The basic properties of the discussed distributions will be collected in the fol- 
lowing Lemma that contains facts from mostly [7] and the paper [3]. 

Lemma 1. i) For n,m > : 

^^^^HAx\q)H„Ax\q)fN{x\q)dx=l^ itZ' 
ii) For n > : 

Hn {x\q) fcN {x\y, p, q) dx = /?"i7„ {y\q) . 

Siq) 

Hi) For n,m > : 

^^^^Pn{x\y,p,q)Pm{x\y,p,q)fcNix\y,p,q)dx^l^ itZ ' 

iv) 

fcN {x\y, Pi,q) fcN {y\z, P2,q) dy = fcN {x\z, p^p^, q) . 

Sig) 



8 



PAWEL J. SZABLOWSKI 



v) For \t\, \q\ < 1 



convergence is absolute, where Wi (q) is defined by \2.()]) . 

vi) For (1 - q)x^ < 2 and V(l - q)t^ < 1 : 

OO 

V{x,t\q) = n {l-{l-q)^tq'' + {l-q) iV')"' > 
fe=0 

convergence V2. 9\) is absolute in t & x and uniform in x. Moreover Lp{x,t\q) is 
positive and /g^^) ^ (x, t\q) fpf {x\q) dx — 1. ip {t, x\l) — exp (^xt — t^/2) . 

vii) For (1 - q) niax(a;^, y^) < 2, \p\ < 1 and V(l - < 1 : 

(i-a-g)p2/^g" + (i-g)pW") 

Tix,t|y,p,gj - JJL {l-{l-q)xtq'^ + {l-q)t^q^'^) ' 

convergence h2.10\) is absolute int&x and uniform in x. Moreover t (x, t\0^ p, <q) is 
positive and J^^^^t {x,t\y, p, q) fcN (x\y.,p,q) dx ^ 1. t {x,t\y,p, 1) = exp (t {x - py) - t^{l - p^)/2) 
via) For (1 — g) niax(x^, 2/^) < 2, |p| < 1 : 

(2.11) fcN {x\y, p, q) = /at ^ -— iJ„(a;|(7)iJ„(2/|g) 

and convergence is absolute t, y & x and uniform in x and y. 

Proof, i) It is formula 13.1.11 of [7\ with obvious modification for polynomials 
Hn instead of /i„ (compare (|2.3p ) and normalized weight function (i.e. /at) ii) 
Exercise 15.7 of [7] also in [1], iii) Formula 15.1.5 of [7] with obvious modification 

for polynomials P„ instead of Pn {x\y, p,q) = (1 — q)"^/^P ^-^^\--^=^ p^q^ and 

normalized weight function (i.e. few), iv) see (2.6) of [3]. v) Exercise 12.2(b) and 
12.2(c) of [7|- vi)-viii) The exact formulae are known and are given in e.g. [7] 
(Thm. 13.1.1, 13.1.6) and |10| (3.6, 3.10). Absolute convergence of iy9 and r follow 
(12. 5p and v). Positivity of ip and t follow formulae 1 — (1 — g) xtq'' + (1 — q) t^q^^ 
= (1 - q){tq^ - x/2Y + 1 - (1 - q)x'^/'i and 1 - (1 - g) pytq^ + (1 - g) p^t^q^^ = 
(1 - q)p'^(q^t - v/{2p)f + 1 - (1 - g)yV4- Values of integrals follow and 
(|2.10p and the fact that {Hn} and {Pn} are orthogonal bases in spaces L [q) and 
CL{y,p,q). □ 

Corollary 1. Every marginal distribution of multidimensional q— Normal distri- 
bution Nd {va.,<y'^ T p\q) is multidimensional q— Normal. In particular every one- 
dimensional distribution is q— Normal. More precisely i~th coordinate of Nd (m, cr^, pjg) — 
vector has {mi,af,q)— Normal distribution. 

Proof. By considering transformation {Xi, . . . , Xd) — > ( ^V™^ ' • ■ • ' ^VJ"'' ^ 
we reduce considerations to the case Nd{0,l, p\q). First let us consider d — 1 
dimensional marginal distributions. The assertion of Corollary is obviously true 
since we have assertion iv) of the Lemma [1] We can repeat this reasoning and 
deduce that all d — 2, d — 3, . . . , 2 dimensional distributions are multidimensional 
g— Normal. The fact that 1— dimensional marginal distributions are g— normal 
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follows the fact that fcN {y\x,p,q) is a one-dimensional density and integrates to 
1. □ 

Corollary 2. //X= (Xi, . . . , X^) - (m, 1, then 
i) yn e N,l < ji < j2 ■ ■ ■ < jm < i < d : 

Xi\Xj^^ , . . . , Xj^ ^ fcN {x,\xj^ , rife^j^ P/c : • Thus in particular 

E {H„ {X, - rrn) \X,, , . . . , X,„J - f ff p}\ Hn {Xj 

and var {Xi\Xj-^ , • • ■ , Xj 
a) Vn € N, 1 < ji < 




Xi\Xj^,. ..Xj^^,Xj^,. . ■,Xj^ Jn {x.i\q) Y[ 



^J-^ hi {Xi ,Xj^,pl,q)hi{xi, Xj^ ,p*m,q)' 



where hi {x,y,p,q) = {{1- p'^q'^^)^ - {l-q)pq'' (l + p^gZ/j xy+{l-q)p'^q^''{x^+y^)), 

P*k ~ WiJj^ Pi ■> Pm — YVi=i ^ Pi ■ Thus in particular this density depends only on 
Xj^, and Xj^ . 

Proof, i) As before, by suitable change of variables we can work with dis- 
tribution Nd (0, 1, p\q) . Then following assertion iii) of the Lemma [Hand the fact 
that m— dimensional marginal, with respect to which we have to integrate is also 
multidimensional g— Normal and that the last factor in the product representing 

density of this distribution is fcN {xi\xj^^ , Jlfc^j^ Pk^ we get i). 

ii) First of all notice that joint distribution of {Xj-^ , . . . Xj^ , Xi,Xj^ , • ■ • , Xj^^) 
depends only on Xj^,Xi,Xj^^ since sequence Xi, i = 1 . . . ,n is Markov. It is also 
obvious that the density of this distribution exist and can be found as a ratio of joint 
distribution of {Xj^,, Xi, Xj^) divided by the joint density of . Keeping 

in mind that Xj^ , Xi, Xj^ have the same marginal /at and because of assertion iv 
of Lemma [T] we get the postulated form. □ 

Having Lemma [1] we can present Proposition concerning mutual relationship 
between spaces L (q) and CL [y, p, q) defined at the beginning of previous section. 

Proposition 2. S (—1,1),?/ G S [q) ,\p\ < 1 : L (q) = CL{y,p,q). Besides 
3Ci iy,p,q) , C2 iy,p,q) : blL < Ci \\g\\cL "."^d \\g\\cL < ^2 UWl for every g £ 
L{q). 

Proof. Firstly observe that : (1 — p'^q^^Y ^ (1 ^ l)pi''{^ + p^q^^)xy -I- (1 — 

2_ ■ 

2 

which is elementary to prove. We will use modification of the formula (|2.1ip that 
is obtained from it by dividing both sides by /jv {x\q) . That is formula: 



q)p^{x^+y^)q^^' = {\-q)p'q^'' [x - y P'^ + (1 - (1 - g)yV4)(l - p2g2fe)2 



Ha 



(1 — p^q^^Y — (1 — q)pq^{l + p^q^^)xy + (1 — q)p^{x^ + y^)q^ 

oo 

Y^^^Hn{x\q)H,,{y\q). 

n=0 ^ k- 
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Now we use p.Sp and assertion v) of Lemma [T] and get Vx, y £ S (q) : 

fcN {x\y,p,q)<.fN{x\q)^^. 

(P)oo 

Hence C2 = and Ij^H^^ < , for every ge L{q). Thus g G CL{y,p,q). 

Conversely to take a function g e CL {y, p, q). We have 

00 > / \g{x)\^ fcN {x\y,p,q)dx. 
Js{q) 

Now we keeping in mind that {l-p'^q^^)'^-{l-q)pq^{l + p^q^^)xy+{l-q)p^{x^ + 
y'^^q'^k quadratic function in x, we deduce that it reaches its maximum for 
X & S (q) on the end points of S (q) . Hence we have 

(1 - pVr - (1 - q)pq'(l + pV')xy + (1 - q)p'{x' + y')q'' 

< (i + pV + v^lralkl'T- 

Since for Vy G S* (q) , \p\ , |q| < 1 : 

00 

+ pV + VT~q\py\\qh'<^ 

and we see that 

00 > \g{x)f fcN ix\y,p,q)dx ^ / \g {x)f /n {x\q) 

JS{q) JS(q} 

xfi ,. ^ 



k=0 



1 - p^q*^ 

(1 — p'^q^^Y — (1 — q)pq^{l + p'^q'^'^)xy + (1 — q)p'^(x^ + y'^)q'^^ 



> J— I \g {x)\^ Jn {x\q) dx. 

So g€ L{q). □ 

Remark 2. Notice that the assertion of Proposition^^ is not true for q = I since 
then the respective densities are N (0, 1) and N (^py, 1 — p^) . 

Remark 3. Using assertion of Proposition\^we can rephrase Corollary\^in terms 
of contraction TZ {p, q) , (defined by \2.12\l . below). For g £ L (q) we have 

E {g (X.) \X,, ,...,X,J=Tz(^ n 9 j (3 (^^■". )) ' 

where TZ {p, q) is a contraction on the space L (q) defined by the formula ( using 
polynomials _ff„ for \p\ , \q\ < 1) ; 

00 00 

(2.12) L{q)3f^Yl ^ ^ (P' 1^ (/) = E ""^P^^^ (^1'^) • 

i=0 i=0 

By the way it is known that TZ is not only contraction but also ultra contraction i.e. 
mapping L2 on Loo (Bozejko). 

We have also the following almost obvious observation that follows, in fact, 
from assertion iii) of the Lemma [1] 
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Propositions. Suppose that ^={Xi, Xd)^ {0,1, p\q) andg£L{q). As- 
sume that for some n G N, and 1 < ji < j2 ■ • • < jm < i < d. 

i) If E {g (Xi) \Xj-^, . . . , Xj^^) —polynomial of degree at most n of Xj^, then 
function g must be also a polynomial of degree at most n. 

ii) If additionally Eg {Xi) = 

(Generalized Gebelein's inequality) E((E(.g(XO|Xji , . . . , X^^)'^) < r'^Eg'^{Xi), 
where r = nl=j„, Pk- 

Proof, i) The fact that E {g (Xi) \Xj^ , . . . , Xj^) is a function of Xj^^ only, is 
obvious. Since g & L{q) we can expand it in the series g {x) — X)i>o ^i^i {x\q). By 
Corollary [H we know that E {g {Xi) \Xj^ X^Jj = Y.^■^^ c,r^H~{Xj^ \q) for r = 
Y[k=j Pk- Now since Cir' = for i > n and r ^ we deduce that q = for i > n. 

ii) Suppose g{x)^J2T=l9^H^{x)■ We have E(5(X,)|Xji , . . . , Xj„J = E^i 5»7-'^^(lS,„). 
Hence E((E(5(X,)|X,,, . . .,X,J^) = ^.^^i 9^ W,! < EZigf - r^^gH^.). 

a 

Remark 4. As it follows from the above mentioned definition, the multidimensional 
q— Normal distribution is not a true generalization of n~ dimensional Normal law 
N„ (m, S). It a generalization of distribution N„ (m, S) with very specific matrix 
S namely with entries equal to = af; aij — ai<7j Ilfe=i Pk f'^''~ * < J '^'^'^ '^ij = 
(jji for i > j where Ui ; i = 1, . . . ,n are some positive numbers and < 1, i = 
l,...,n~l. 

Proof. Follows the fact that two dimensional g— Normal distribution of say 
{Xi/(j.„Xj/aj) has density f^ {xi\q) fcN (a^iki, lid Pk^q) if « < J- ^ 

Remark 5. Suppose that {Xi , . . . , Xn) ^ Nn (m, cr, p\q) then Xi , . . . , Xn form a fi- 
nite Markov chain with Xi ^ {mi, af, q)— Normal and transition density Xi\Xi-i — 
y ^ {rui, crf,y, Pi^i, q)- Conditional Normal distribution 

Following assertions vi) and vii) of Lemma [1] we deduce that for Vi^ < 1 / ( 1 - g) , 
IpI < 1 functions {x,t\q) fi^ {x\q) and t {x,t\y, p,q) fcN {x\y, p,q) are densities. 
Hence we obtain new densities with additional parameter t. This observation leads 
to the following definitions: 

Definition 2. Let \q\ G (— < 1/(1 — q),x G S{q). A distribution with the 
density if {x, t\q) f^ {x\q) will be called modified (t. q) —Normal (briefly {t, q) ~MN 
distribution). 

We have immediate observation that follows from assertion vi) of Lemma [TJ 

Proposition 4. i) /^^^^^ (p {y, t\q) fN {y\q) fcN {x\y, p, q) dy ^ ip {x, tp\q) fN{x\q) 
ii) Let X - {t, q) -MN. Then /or n G N : E (iJ„ {X\q)) = t". 

Proof, i) Using assertions vi) and viii) of the Lemma [1] we get: 
/ y^{y,t\q) fN {y\q) fcN {x\y,p,q)dy 

JS{q) 

= fN{x\q) fN{y\q)y2TTi^'(y\'i)'l2TTi^i^y\'i^^j^^\'i^^y- 
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Now utilizing assertion ii) of the same Lemma we get: 

^ iy, tW) In {y\q) Icn {x\y, p, q) dy 



S{q) 



oo 



= In {x\q) ^ -TTtH^ {x\q) = LP (x, tp\q) fN{x\q). 

To get ii) we utilize assertion vi) of the Lemma [1] □ 
In particular we have: 

Corollary 3. If X ^ {t,q) - MN, then EX = t, var (X) = 1, E (^{X - tf^ = 
-t{l^q),E(X-t)^ = 2 + q-e{5 + 6q + q^) . 

Proof. We have a;^ = Hi, {x\q) + {2 + q) Hi {x) and H^ {x\q) + (3 + 2(7 + 
q^)H2 {x\q) + 2 + q = x^ so E{X - tf = E {X^) - 3E {X^) t + 3E {X) ^ - 
= t^ + {2 + q)t-M{l + t^) + 3t3 -t^ = -t{l~ q) and E{X-tf = E (X^) - 
4tE {X^)+QeE (X2)-4t3E = t^-{Z+2q+q^)e+2+q-At (t^ + {2 + q)t) + 

(<2 + l) - At"^ + which reduces to 2 + q - + Qq + q^) . □ 

In particular kurtosis of — MN distributions is equal to — (1 — g) — 

(5 + 6g + (7^). Hence it is negative and less, for t 7^ 0, than that of g— Normal 
which is also negative (equal to — (1 — 5)). 

Assertion i) of the Proposition |3] leads to the generalization of the multidi- 
mensional (7— Normal distribution that allows different one-dimensional and other 
marginals. 

Definition 3. A distribution in having density equal to 

V {xi,t\q) Jn {xi \q) Y[ fcN {xi+i \xi,p^, q) , 

i=\ 

where X, e S (q) , Pi e (-1,1)\{0}, i = 1, . . . , d-1, \q\ e (-1,1], < l/(l-g) will 
be called modified multidimensional q-Normal distribution (briefly MMNd{p\q, t) ). 

Reasoning in the similar way as in the proof of Corollary [T] and utilizing obser- 
vation following from Proposition^ we have immediately the following observation. 

Proposition 5. Let {Xi, . . . ,Xd) ^ AIMNd{p\q,t). Then every marginal of it 
is also modified multidimensional q— Normal. In particular yi = l,...,d : Xi ^ 

{tlYk=\Pk,Q)-MN 

Remark 6. Suppose that {Xi, . . . Xd) ^ MMNd{p\q, t) and define Pq = 1, then the 
sequence Xi, . . . Xd form a non- stationary Markov chain such that Xi ^ Lp (^\t Ofe=i Pk^ In 
with transitional probability XiJ^i\Xi ~ y ^ fcN {x\yTPi,q) ■ 

We can define another one-dimensional distribution depending on 4 parameters. 
We have: 

Definition 4. Let \q\ £ (—1, 1], < 1/(1 — q),x,y £ S (q) , \p\ < 1. A distribution 
with the density r (x, t\y, p, q) fcN ix\y, p, q) will be called modified {y, p, t, q) —Conditional 
Normal (briefly {y, p, t, q) -MCN). 
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We have immediate observation that follows from assertion vii) of Lemma [TJ 
Proposition 6. Let X - (y, p, t, q) -MCN. Then /or n G N : E (P„ {X\y, p, q)) = 

Hence in particular one can state the following Corollary. 

Corollary 4. EX = py + {\ - p'^)t , var(X) = {\- p'^)[\-{\- q)typ+{\- q)t^ p'^). 

Proof. Follows expressions for first two Al-Salam-Chihara polynomials. Namely 
we have: Px{x\y,p,q) = x- py, P2{x\y,p,q) = x"^ -1 + p^ + qp'^y'^ -xpy[l + q). □ 

We can define two formulae for densities of multidimensional distributions in 
R'^. Namely one of them would have density of the form 

'^{xi,t\q) /at {xi\q)T{x2,t\xi,p^,q) J| JcN {x,+l\x^, p^,q) 

i=l 

and the other of the form 

d-l 

/at {xi\q) T (x2,t\xi,pj^,q) Y[ fcN {xt+i \xi,p^, q) . 

1=1 

However to find marginals of such families of distributions is a challenge and an 
open question. In particular are they also of modified conditional normal type? 

3. Main Results 

In this section we are going to study properties of 3 dimensional case of mul- 
tidimensional normal distribution. To simplify notation we will consider vector 
(y, X, Z) having distribution N3((0, 0, 0), (1, 1, 1), (Pi, /52)k) that is having density 
fcN {y\x,PiTq) fcN{x\z, p2,q),fN{z\q)- We start with the following obvious result: 

Remark 7. Conditional distribution X\Y,Z has density 

fo,\ j,f I \ f I \ ^(^l'^2)ooTT Wk{y,z,p^p2,q) 

(3.1) (I){x\y,z,p^,p2,q) ^ fN{x\q) , 2 2\ \\ 7 ^ 7 T' 

(PiP2)oo r=o i^^y^Pi^1)Wk [x,z,p2,q) 

where we denoted Wk {s,t,p,q) = (1 — p^q^'^Y - (1 — q)pq^{\ + p'^q^^)st + (1 — 
q)p^{s^ +t^)q^^. 

Proof. It is in fact rewritten version of the proof of assertion ii) of Corollary 

m □ 



Remark 8. Notice that (j) {x\y, z, pi, p2, q) is the re-scaled Askey-Wilson density. 

2 



Namely cj) {x\y^ z, Pi,P2,q) — V'l 2^^!'^' ^' ''^^^^^ 



a = -^p,{y-^^J^^-y^),b^J^p,iy + ^^^^-y^), 



^ ^-Piiz - i\l -r- = '^ Piiz + i\l z'^). 



2 ' \l 1-q ' 2 V 

and ^{t\a, b, c, d) is a normalized (that is multiplied by a constant so that its integral 
is 1) weight function of Askey-Wilson polynomials. Compare e.g. [SJandj?]. Hence 
our results would concern properties of Askey - Wilson density and Askey- Wilson 
polynomials. 
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Let us denote 

9niy,z,p^,P2,q) ^ / Hn{x\q)(j){x\y,z,p^,p2,q) . 

Js{q) 

Our main result is the following 

Theorem 1. i) (j) {x\y, z, p^^, p^lq) = Jn {x\q) x 

Yln=o "[lif 9n iy,z,pi,p2,q) , whcre g„ is a polynomial of order n in {y,z). 

a) More over polynomial <;„ has the following structureyn > I : gn (jj, z, pi, p2, q) 
= Er=o ["] qPiP2^'^t.ri-i {y,z,pj^p2\q) , where Qkd {y,z,p-^p2\q)isa polynomial in y 
of order k and in z of order I. Moreover Oq « [y, z,0\q) ~ Hn iz\q) and 0„ o [y, z, 0, q) 
= Hn {y\q) ■ 

Remark 9. Assertion i) of the Theorem]^ is in fact a generalization of Poisson- 
Mehler formula (that is assertion viii) of Lemma\^ for Askey-Wilson density. 

Remark 10. Notice also that for q — 1, (p is a density function of normal distri- 
bution N I — — — - — J— ^ I (ind it is obvious that expectation of 
\ ^ P1P2 ^ P1P2 j 

any polynomial is a polynomial in — — ''i^ — • Hence it turns out that this 

Pi P2 

is true for all q^ Normal distributions for q G (—1, !]• 

As a Corollary we have the following result. 

Corollary 5. Let X= {Xi, . . . , X^) ^ (0, 1, p\q)- Let us select indices 1 < ji < 
■■■jk <i < jrn < ■ ■ ■ < jh < d. Then 

(3.2) Vn e N:E{Hn{X,\q)\X,„...,X,„X,^,...,X,,)^ 

ln/2\ n-2r 

E E 4tln/2i+r+lHl{X,Jq)Hn-2r-l{X,Jq), 
r=0 1=0 

for [^^J [^^2^J constants (depending only on n, q,p and numbers i, jk jm) 
r = 0, . . . , [n/2\ , s = - [n/2\ + r, . . . , - [n/2\ +r + n-2r. 

Corollary bellow gives detailed form of coefficients aI"] for n = 1, . . . , 4. 

Corollary 6. Let X= {Xi,..., Xj.) (0, 1, p\q). Let 1 < i - 1< i < i + 1 < d. 

Then: 

E(i/„ (X,) = ELT4:-Ln/2j+r+;^' iX,^^\q)Hn-2r-l (X,+i|g) , 

where ^o,-L"/2j+; = lilg (pCipI) ' I ^ 0, . . . ,n, n = 1 . . . ,4. If n 

< 3 then 4"-L«/2j+/ = " ~ ^Ig /'«-i/'«4"-Ln/2j+/' ^ ^ ■ ■ ■ ,n - 1, If n ^ 4 
then aQ = -m,P^-iP.A^£, J = -1,1 and ^^^^ ^ -m^P.^^pAt 4'^ = 
q{l + q)pt,p^Ai% 

In particular : 

il-pf_,){l-pi) 



(3.3) var(X,|Xi,...,X,_i,X,+i,...,Xrf) 



(3 4) y. (1 ~ ^^'^^ ~ P^-lP^X^+l) {^i ^ -^»-lP^-lP^) ^ 

(1 - PtlPI) 
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Remark 11. Notice that in general conditional variance 

var . . . , Xi-i, Xi+i, . . . , Xd) is not nonrandom indicating that q— Normal 

distribution does not behave as Normal in this case, however if we set q — \ in 

Ii3.3\) then we get var . . . , Xi^i, X^+i, . . . , Xd) ~ — / '"^^ \ — as it should 

be in the normal case. 

Notice that examining the form of coefficients for n = 1, ... ,4 we can 

formulate the foUowing Hypothesis concerning general form of them: 

Conjecture 1. For n > 1, we have: Af,"-* , /oi , , = [?] s — ' , I = 

•> — ' 0,-Lri/2J+/ lUq (pi_iPi) 

0, . . . , ri. Moreover for r = I n/2 I and I ~ 0, . . . , n—2r A*'"'' , , , , MI"'' i mi , , 

= Pi^iPiQr,i (<z) , where Qr.i {q) is a polynomial in q with coefficients depending only 
on r and I. 



4. Proofs 

Proof of the Theorem [1] is based on the properties of the following function 
Gi,k (y, z, t\q) = Em>o T^.^rn+i (yl?) ^m+fe {z\<l) ■ Wc wiU nccd somc of its prop- 
erties. Namely we will prove the following Proposition which is in fact a general- 
ization and reformulation (in terms of polynomials i?„) of an old result of Carlitz. 
Original result of Carlitz concerned polynomials Wn{x\q) — X]"=o ['^Iq-''' ^'^'^ 
pressions of the form X^i^o '"''■^1'^-'™'+''*-'"!^-** , compare [3, Exercise 12.3(d) or [6]. 

Proposition 7. i)\/k,l>Q : Gk,i {y, z, t\q) = Gi^k {z, y, t\q) 
ii) for I < j < k : 



(4.1) Gu{y,z,t\q) - J2^-iy 



qi^fHk-^iy\q)Go^;+l{y,z,t\q) 



-i-iyqii)Y^ ^ t^Gk-^,^+l{y,^^,tk)■ 



Hi) 



i=0 



q(^^)fHk-^{y\q)Go4y,z,t\q). 



(4.2) Gk,oiy,z,t\q) = J2i-^y 
iv) 

(4.3) Gk.o {y, z, t\q) = _ J 

^ 1=0 

whereTk,i {y,z,t\q) {Hk-i{y\q)Go4y, z,t\q)+{-l)''q&t''Hk-i{z\q)Gi,o{y, z,t\q)). 

Proof, i) is obvious. 

iii) Take j — k and ^ = in ii). 



fc-i 



t'Tk,t {y,z,t\q) , 
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ii) To prove (|4.ip we will use formula 
iJ„+™ {x\q) = iJ„ {x\q) H^{x\q) - J^fil^'^'"'^ [;] ["] [j],!i/„+™_2,(x|(z). We have 



Gk,i [y, z,t)^ V —^H^+k {y\q) H^+i {z\q) 

^ — ' TO ! 

m>0 ^ ^1 

[i]q^-Hk+m~2j{.y\q))Hm+l {z\q) 

Hk{v\q)Go^i{y,z,t)~^t^ X! i (2^lg)-^'+J+("'-j) 

Gfc-j,i+j(?/,z,t). 



inin(fc.m) 

E 



'k 




TO 


.i. 


9 


.i. 



i=i 



Hence let us assume that ()4.ip is true for j = 1, 2, . . . , to. We have after applying 
just obtained formula (for Gk,i) applied for k~ > k — m and I— > m + 1 : 



m — l 



Gk,i{y,z,t) = Y.^-iy 



1=0 



g(^)fi?fe_,(y|g)Go,.+,(y,z,t) 









i 


q 


TO — 1 



fGk-i,i+i{y,z,t) 



i=0 



g(^)fi/fc_,(y|g)Go,.+,(y,z,t) 



+ 



t'^{Hk-miy\q)Go^m+i{y, z, t) 



fe-r 



E 



k — m 
i 



fGk — m—i.l-\-m-\-i 

{y,z,t)) 







■ i - 1 " 




9 


TO — 1 



fGk-iA+i{y,z,t) 
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Nowsince['=] -H ['-^1 ([M -['-^1 ^^1 and (™) 



-\- VI — {"^2'^) we have 

m 

Gfc,z(y,z,i) = ^(-1) 



1=0 



g(2)i^i7fc_,(y|g)G'o,i+/(y,^,i) 



"fc' 




k — m 




'k 




' i-l' 


m 


9 


i — m 


1 


i 




TO — 1 



)fGk-^,i+t{y,z,t) 



q(^^h'Hk-.^iy\q)Go,^+l{y, z,t) 



i—m-\-l 



'k' 




'i-l 


i 


q 


m 



iv) For fc = this is obviously true. Now let us iterate (|4.2p once, applied however, 
for Go^fc. We will get then 

Gfc,o (y, z, t\q) = E -=0 (-1)' ['] ,9(^)i^^fe-.(yk)Go.(2/, 2, g^'^^'^Go,,. {y, z, t\q) 

(-1)' ^('^^'^(E to (-1)' [ ■] ,g(^)i^i?fc-.(^k)G.,o(y, z, i|g))+9'=('=-i)i2^-Gfe.o (y, z, t\q). 
Thus we see that since for all i < fc — 1 Gi^ and Go,i and are of the claimed form 
then from (14. 3p it follows that Gk o has the claimed form. □ 



Now we are ready to present the proof if Theorem [TJ 



Proof of the Theorem [T] To prove i) we will use formula viii) of Lemma 
[TJ that is Poisson-Mehler expansion formula. Following p. II) we see that 

oo ^ oo ^ 

<?!)(x|y,z,pi,p2,g) = fN{x\q)Y,-fTlH^{x\q)H,{y\q)xY,-fr]H^{x\q)H,{z\q) 

1=0 i=o '-^'q' 

oo ^ ^ 

z=o ^ ^q- 

First, let us concentrate on the quantity: 

R{x,y,z,p^,P2\q) = TT-jHn ix\q)Hn iy\q) x ^ fTl^^ ix\q) {z\q) . 

n>0 ^"^9- m>0 ^'^^Q' 

We will apply identity (12.41) . distinguish two cases n + to is even and n + to is odd, 
denote n + to — 2j — 2fc or n + to — 2j = 2fc + 1 depending om the case and sum over 
the set of {(n, m) : n + m — 2k < 2 min(n, to), m,n> 0}U{(n, to) : n + to — 2A; — 1 
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< 2 min(n, m), m, n > 0}. We have 
R{x,V,z,p-^,P2\q) = ^ . /, . , gn {y\q)Hrn {z\q) 
H2kix\q) 



n 




m 


J. 


9 





E 



■EE 



E 



^ [2fcl ! ^ ^ [m - 7l ! \2k - (to - j)! ! hi 

„ . I ^ oo 2fe+l+j rol, I 11 I 2fe+l-i+j 



E E 



j] '[2fc + l-(m-j)] ' 



i7„i (yl g)iJ2fe+2:,-m 1 g) 

i/m (y\q) H2k+i+2j-m {z\q) 



E 



^^2fc(2:|g) 



/c=0 ^ ^9 i=0 1=0 ^-^^9 L J9 L J9 



^ (:r |g) ^ [2k + 



k=0 ^ '1 j=0 i=0 

We get then 



{y\q)H2k+i+j-t {z\q) 



2k 



R{x,y,z,p„P2\q)^}_^ 2^ 

fe=0 ^^'^^9- i=0 



P1P2 2^ 



■E 

fe=0 



H2k+i{x\q) 
I ! 

'9 j=0 



2/c+l 

E 



2fc + 1 



i 2k- 
P1P2 



(P1P2P 



g« {x\q) -A 

„=0 ^ ^9- 1=0 



j=o ^■'^9 
{PiP2y 



^1^2 2^ r . , 
J j=0 ^■^^9' 



H^+J{y\q)Hr,-^+j{z\q) 



Usmg mtroduccd in Proposition[7]function d^k [v, z, t\q) = J2m.>o |^-^™+/ ivW) Hm+k [z 
we can express both 

E^o {y\q)H, {z\q) = Go,o{y, z, PiP2\q) and 



R{x,y,z,p^,P2\q) = E„=. 



r=o^T^ELo [:]^p\P2-'G.,n-^{y,z,p,p2\q). Our The- 
orem will be proved if we will be able to show that V/,fc > : Gi^k {y, z,t\q) = 
Go,o [y, z, t\q) Qi^k {y, z, t\q) where Qi^k is a polynomial of order a Hn y A: in 2. This 
fact follows by induction from formula ()4.3p of assertion iv) of the Proposition [7] 
since it expresses Gkfi in terms of k functions G^^o and Go,i for i = 0, . . . , fc — 1 and 
the fact that all Gi^k can be expressed by Gi^ and Go,i ; i < k + I. □ 

Proof of Corollary [5l By Theorem [1] we know that regression 
E(iJ„ {Xi\q) \Xj-^^ , ■ ■ ■ , Xjf_ , Xj^ , • ■ • , Xj^ ) is a polynomial in Xj^ and Xj^ of order 
at most n. To analyze the structure of this polynomial let us present it in the form 
Sr=o '^s,n-ffs {Xj^\q) where coefficients as^„ are some polynomials of Xj,^. Now let 
us take conditional expectation with respect to Xj-^ , . . . , Xj^ of both sides. On one 
hand we get 



E{H„iX^\q)\X, 



.,x, 



\''n=3k 
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on the other we get 

/ i^-i \ ' 

^fls.n H ^" Hs{Xjk\q). 
s=0 \m=jk / 

Since as,n are polynomials in Xj^, of order at most n, we can present them in the 
form 

n 
t=0 

Thus we have equality: 

i-l \ " " / ^""^ 

Now we use the identity (I2.4p and get 

n I Un (^j.k) = EE/^M I n 



inin(t,s) 

5 



E 



X 

m=0 



m 



9 
9 



Hence we deduce that (3^ ^ — for t + s > n , t + s = n — l,n — 3, . . . To count 

the number of coefficients A^'^ observe that we have n + 1 coefficients Ag"^ since 

k ranges from — [^J to — [^J + n , n — 1 coefficients A^"^ where k ranges from 
- [n/2\ + 1 to - [n/2j + n - 1 and so on. □ 

Proof of Corollary [H The proof is based on the idea of writing down 
system of [^^J [^^J ( n — 1,...,4 ) linear equations satisfied by coefficients 

fe- These equations are obtained according to the similar pattern. Namely 
we multiply both sides of identity p.2p by (Xi^i) and Hk (Xi) and calculate 
conditional expectation of both sides with respect to J-^i or with respect to J^>i 
remembering that E = pT^^pTHrr. (^^-l) and E (i?™ {X,) 

= pl"iiHm (Xi^i) and similar formulae for J-yi. We expand both sides with re- 
spect to Hs, s = n + m + k — 2t,t = 0,..., [{n + m + k)/2\ . On the way we utilize 
and compare coefficients standing by on both sides. Thus each obtained 



equation involving coefficients A^^^j can be indexed by s,m,j and r if we calculate 
conditional expectation with respect to J-^i of I if we conditional expectation is 
calculated with respect to J-yi- Of course if s = then r and I lead to the same 
result. Formulae for Aj-"'', for n = 1 are obtained by taking s = 1, m — 0, j = 
and applying r and /. For n — 2 first we consider m — 0, j — and s = 2 
and applying r and I and then m = 0, j = and s — 0. In this way we get 3 
equations. The forth one is obtained by taking m — 0, j — 1, s — 1 and r. Denote 
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X = {aI^_j^, Af^l, A''qI,A'{^q)'^, then X satisfies system of linear equation with, ma- 



l(2) a{2) .(2) 



( 



trix 



1 

pi-iPi 





pi-iPt 

Pt~lPz 
Pt-lPz 



Pi-lpi 
1 





( 



and with right side vector equal to 



pU 
pI 






1 

\ 



Besides formulae for coeffi- 



cients A^fc, for n 





V \\P^-^P^ 

1, 2 can be obtained from formulae scattered in the literature 



for n 



like e.g. [T], |12j or jl9j . To get equations satisfied by coefficients A^ .^ 
3, 4 First n + 1 equations are obtained by taking m = , fc = and s =3, 1 if n = 3 
and s =4, 2, if n = 4 and then applying operations r and I. Then, in order to get 
remaining 2 (in case of n = 3) or 4 (in case of n = 4) equations one has to be more 
careful since it often turns out that many equations obtained for some m and fc are 
linearly dependent on the previously obtained equations. In the case of n = 3 to 
get remaining two linearly independent equations we took m = 2, fc = 0, s = 3 and 
applied operations r and /. In this way we obtained system of 6 linear equations with 
/I 

Pi-iP 

(1 

{l + q)ptiPj 

i + mlpiipi 
m,p^-lP^ + mgPtlPl i + [2]>f_ip: 



matrix 



[3; 



pI-ipI 




Pt 



ipt 

V maPtiPi 



2 2 

Pi~iPi 

i^+q)Pt^lPr 
['^],P^-lP^ + [3']qP|-lP! 



Plipf 

1 

1 

[^]nPt-^lPi pI-iPI 



( 



right hand side vector 



.(2) .(3) J(3)^T 
■ ^0,2' ^1,0' ■ 



pU 
pI 




[^]qPl-lP^ 
V [^]gP^-lP^ J 



if the vector of unknowns is the following 



For n 



4 remaining 4 equations we obtained by 
taking: (m = 1 , k — 4, s— 3, r), (m = 4, fc = 1, s= 1, r), {m = 4, fc = 2, 
s = 4, r) and (m = 2, fc = 4, s = 4, , r). Recall that in this case we have 9 
equations. Matrix of this system has 81 entries. That is why we will skip writ- 
ing down the whole system of equations. To get the scent of how complicated 
these equations are we will present one equation. For n — 4 one of the equa- 
tions (referring to the case m = 4, = 1, s = 1, r) is [2]^ (l + q^) [3]^ ([4]^ -I- 

[5],))ft-ift4!-2 + [2]^ (1 + m,x{l + ptiPl + SqpUpH3q^PtipHQ^ptiPl 



+ q^pl_^pl)Al 

^3 „2 



(4) 



+ [2]^(l + g2) {Z\p.,_^pm, + pUpI 



q-pU.pl + q^pUpDA'-^l + [2]^ (1 + q^) [3], p\^p2^{Xi\ + ([4]^ + [5]^)pt ^pf + 



[21,(1 + 9^) [3],plV2^(l 



\pt 



q^pI-ipI 



<fpLiPi)A^- 



[2]^(l + 'z')[3],Pti([4], + [5],p|_i)p,. □ 






Pi-iPi 
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2 2 

pf-iPf 
Pi-iPi J 
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